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The role of the operator-product expansion in QCD calculations is discussed. Approximating the two-point 
correlation function by several terms and assuming an upper bound on the truncation error along the euclidean 
ray, we consider two model situations to examine how the bound develops with increasing deflection from the 
euclidean ray towards the cut. We obtain explicit bounds on the truncation error and show that they worsen with 
the increasing deflection. This result does not support the conventional believe that the remainder is constant for 
all angles in the complex energy plane. Further refinements of the formalism are discussed. 



1. Introduction 

A standard way to investigate the interconnec- 
tion between perturbative and non-perturbative 
effects in quantum chromodynamics is the 
operator-product expansion (OPE) M. In the 
momentum space, it has the form (we put y = 0) 

i [ dx^ x A(x)B(0) a J2 C k(q)O k , (1) 
fc 

where q is the total four-momentum, q 2 = s = 
—Q 2 . The coefficients Ck{q) are ordered accord- 
ing to the increasing exponent k in s~ fc . 

Unlike perturbative expansions (in powers of 
the coupling constant a s ), very little is known 
about the mathematics of the operator product 
expansion. The knowledge is not sufficient to 
give (0) precise mathematical meaning, not even 
is known the exact composition of terms on the 
right hand side of (|l]) (in particular, whether 
terms exponential in the variable Q^ 1 have to be 
added to describe strong- interaction processes). 
There are reasons to believe that the series is 
divergent Q, but it is not known whether it is 
asymptotic to the function expanded. The large- 
order behaviour of its terms is not known either. 
Last not least, the symbol ~ in (|l]) is understood 
differently in different contexts. 
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In spite of serious mathematical uncertainties, 
the operator product expansion is of paramount 
importance for solving a number of practical 
problems in quantum chromodynamics. To cal- 
culate many measurable quantities in QCD, it is 
necessary to investigate the properties of the OPE 
in the complex Q 2 plane away from the euclidean 
region, i.e., away from the semiaxis Q 2 > 0. This 
leads us, with regard to (|l|), to the problem of 
finding conditions under which a power expan- 
sion of the type 

/(1/Q 2 )«$>(g)/Q 2fc (2) 

can be extended to angles away from the eu- 
clidean semiaxis in the Q 2 plane. But such an 
extension is a delicate problem requiring precise 
mathematical conditions, which are not known 
here. 

In practical applications, it is usually assumed 
that the properties of the operator product ex- 
pansion away from the ray Q 2 > are the same 
as those along it, except for the cut Q 2 < (the 
Minkowski region). Instead of making this sim- 
plifying technical assumption, we prefer to intro- 
duce some model assumptions which would free 
the scheme from a priori neglecting the influence 
of the cut on the truncation error along a ray, keep 
the model possibly close to real situations and, 
simultaneously, give the problem precise mathe- 
matical meaning. 



2. The model 

I will discuss the problem on the example of the 
following model situation. Let F(s) be holomor- 
phic in the complex s plane cut along the pos- 
itive semiaxis, with the possible exception of a 
bounded domain around the origin. Let the num- 
bers a k , k — 0, 1, 2, ...ft— 1, and a positive number 
A n exist such that the inequality 



\F(s)-J2ak/(~s) k \ < A n /\s\ n 



(3) 



fe=0 



holds for all real s smaller than a negative num- 
ber. The problem is under what conditions this 
inequality can be continued to the complex plane, 
and what will be its form away from the original 
region of validity, i.e. away from the negative real 
semiaxis (euclidean region). 

This problem is of interest for QCD because 
it could, when solved, tell us how the inequality 
(||) develops when the ray along which the opera- 
tor product expansion is studied departs from the 
euclidean region towards the Minkowski one. We 
shall discuss the problem under two simplifying 
assumptions. 

Note first that the inequality (g) amounts to as- 
suming that the n-th order remainder R n (—l/s), 



R n (z) = f(z) - } 4 a k z- 



-l 

k=0 



(4) 



tends to zero for z — ► as the n-th power of z for 
at least one value of n. (Here, we use the notation 
z = e lip = —1/s and f(z) = F(s).) This assump- 
tion may appear as too optimistic; indeed, it is 
difficult to make a realiable estimate of the trun- 
cation error even in the euclidean region, because 
it is not known whether there is a limit of QCD 
in which the operator-product expansion becomes 
exact. 

A second assumption we make is that the log- 
arithmic energy dependence of the expansion co- 
efficients dfe can be neglected. 

We do not demand, on the other hand, that the 
series on the left-hand side of (§) be convergent or 
asymptotic to F(s): our approach is more general 
and can be applied whenever the remainder obeys 
(0) at least for one value of n. 



Assuming the bound (||), we observe how it 
varies with the deflection of the ray away from 
euclidean region. It is natural to expect that the 
solution to the problem depends on additional as- 
sumptions imposed on the function F(s) = f(z). 
In particular, the form of the dependence of the 
resulting bound is determined by the assumed 
character of the discontinuity along the cut. 

We consider two possible sets of such additional 
assumptions. 

3. Results 

(i) In the first of them, we assume that f(z) 
admits the integral representation of the form 



p(t) 



zt 



dt 



(5) 



with z complex and p(t) nonnegative for t > 0, 
and that the moment 



„„ = I t n p(t)dt 
lo 



(G) 



exists for at least one positive integer n. Then 
it follows that the remainders Rk(z) as defined 
by (|J), with k = 1,2, ...,n, are bounded by the 
inequalities 



\Rk(z)\ <a k \z\ k 
and 

\R k (z)\<a k \z\ k /\ sin (p\ 



(7) 



(8) 



for Kez > and Rez < respectively. We 
see that the denominator on the right-hand side 
of the inequality (^|) produces a looser estimate 
along rays lying in the halfplane Re z < 0, where 
the cut is present. 

Note that the integral representation (||) is not 
equivalent to the dispersion integral representa- 
tion that is typical for the QCD two-point func- 
tions, because these quantities usually obey a 
once-subtracted dispersion relation, in which case 
the moments (||) do not exist. Contrary to this, 
the approach based on (||) and (||) consists in as- 
suming that at least several lowest moments a n 
are convergent integrals. This feature may be of 
interest in some model situations, which we do 



not consider here. To get closer to physics, we 
consider an alternative scheme below. 

(ii) In this second scheme we assume, instead 
of the above conditions, that f(z) is bounded by 
a constant M inside a circle of radius d in the cut 
z plane. The resulting bound on the remainder 
has the following form: 

\Rn{re^)\<M n (r/d n ) n ^^\ (9) 

Here, M n and d n are constants which related to 
M and d. This inequality is a special consequence 
of a general theorem, which will be published sep- 
arately j|. The inequality ^) gives an upper 
bound on the remainder R n (z) along every ray 
tending to the origin (i.e., to infinite energy), the 
estimate becoming worse with increasing deflec- 
tion from the positive real semiaxis (euclidean re- 
gion), i.e., with the ray approaching closer the 
cut. 

Since the estimate (||) might seem rather loose 
(note that n in the exponent on the right hand 
side is multiplied by an angle-dependent factor 
that vanishes on the cut), it is perhaps worth 
mentioning that there are functions that saturate 
it. Details will be discussed in ||. 

4. Discussion 

The bounds obtained are related to the prob- 
lem of finding error estimates for the QCD calcu- 
lations, in which contour integrals of the type 

67ri I (1 - s/m 2 ) k {\ + 2s/m 2 T ) l P(s)^ (10) 

J\s\— m 2 m 

occur (with k and I being nonnegative integers). 
One example is the determination of a s (m) from 
the t lepton hadronic width ||; in this case, m is 
the r-lepton mass, k—2, 1=1, and P(s) is a com- 
bination of the electromagnetic two-point corre- 
lation functions n(s) for the vector and the axial 
vector colour singlet massless quark currents, the 
II(s) being defined as 

W v = (g^q 2 - ^V)IT(- g 2 ) 




-^(0|T(^(x)r(0))|0) 



for two currents and f . Note that the inte- 
gral ( |ici| ) may receive essential contributions from 



the region near s — m 2 , where the OPE has lit- 
tle chance appropriately to represent the function 
expanded. A fortunate circumstance is that the 
double zero of the kinematic factor (1 — s/m 2 ) 2 
in the integrand suppresses the contribution from 
this dangerous segment. But a quantitative anal- 
ysis of this argument is, to the best of our knowl- 
edge, still lacking. 

It has been the motivation of the present re- 
search to find an explicit quantitative estimate of 
the truncation error. In considering the problem 
of the evaluation of the contour integrals of the 
type (|l(]) , we propose a model scheme of assump- 
tions with the aim to evaluate the influence of the 
cut Q 2 < by estimating the truncation error 
along different rays in the complex Q 2 plane. In 
this scheme, the property (||) occupies the central 
position. Its combination with analyticity allows 
an extension of (|5|) into the complex s plane. 

We have cosidered two possible additional as- 
sumption sets, (i) (||) and (^), and (ii) the as- 
sumption that f(z) is bounded inside a circle of 
the cut z plane. While the case (i) leads to the 
bounds (Q) and (||), the case (ii) yields the bound 
@, in which the exponent itself depends on the 
phase ip, thereby representing a significant wors- 
ening of the asymptotic bound along rays that are 
close to the cut. 

In case (i), an improvement of the integrated 
error in the integrals of the type (|Io|) may be 
reached by reversing the order of the s-integration 
and the z-integration. No such chance seems to 
exist in the case (ii). 

Our results can be considered as a step out of 
the conventional scheme, in which the pragmatic 
apriori assumption is made that the truncation 
error is the same along all rays except the cut. 
Our bounds are based on a plausible, but still 
crude and fictitious picture of the operator prod- 
uct expansion, possibly good for certain model 
situations. In application to physics, however, as 
was pointed out above, the scheme is not fully ap- 
propriate either, due to the simplifying assump- 
tions discussed in Section 2 of this talk. A further 
refinement is desirable. Next step should include 
a regard to the logarithmic g-dependence of the 
coefficients ak(q)- Work along these lines is in 
progress. 
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